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CURRENT  ENHANCEMENT  FOR  HOSE-UNSTABLE  ELECTRON  BEAMS 


1 .  Introduction 

The  propagation  of  an  energetic  charged  particle  beam  is  controlled  by 
the  electromagnetic  fields  and  forces  which  act  on  the  beam.  These  fields  and 
forces  are  modified  by  the  electrical  properties  of  the  medium  surrounding  the 
beam.  An  assessment  of  these  modifications  is  required  to  accurately  predict 
beam  advance.  The  net  current,  defined  as  the  sum  of  the  axial  beam  current 
and  the  axial  conduction  current  Induced  in  the  medium,  is  a  useful  indicator 
of  such  modifications. 

Beams  passing  through  a  gas  ionize  the  gas  and  thereby  change  its 
conduction  properties.  In  low-pressure  gases  a  two-stream  instability  can 
arise  which  propels  plasma  electrons  forward.1  The  forward-moving  plasma 
electrons  raise  the  net  current  for  negatively  charged  beams  but  lower  it  for 
positively  charged  beano.  For  relativistic  electron  beams,  net  currents  as 
large  as  three  times  the  beam  current  have  been  observed. 1-3  These 
enhancements  are  accompanied  by  microwave  emission  characteristic  of  the 
plasma-electron  oscillation  frequency. 

At  higher  pressures,  binary  collisions  between  the  plasma  electrons  and 
the  gas  ions  and  molecules  damp  the  two-stream  instability  and  its  associated 
current  enhancement.  Axial  plasma  currents  then  develop  primarily  from 
inductive  effects  which  typically  drive  plasma  return  current  opposite  to  the 
beam  so  as  to  minimize  changes  in  magnetic  flux.  Met  currents  at  high  gas 
pressures  are  thus  anticipated  to  be  less  than  the  beam  current.1 
Manuscript  approved  February  7, 1986. 


Recent  experiments  performed  on  the  Experimental  Test  Accelerator  (ETA)-3 
and  elsewhere®’^  have  shown,  however,  that  the  net  current  can  exceed  the  beam 
current,  by  as  much  as  a  factor  of  two,  even  at  gas  pressures  of  a  few  hundred 
torr.  The  absence  of  microwave  emission  at  gas  pressures  above  10  torr 
indicates  that  this  enhancement  is  not  due  to  the  two-stream  instability.  A 
second  proposed  explanation  is  energetic  secondary  electrons  (delta  rays) 
swept  forward  by  the  beam  self -magnetic  field.  Estimates  of  delta-ray 
currents  show  that  they  are  too  small  (less  than  20  percent  of  the  beam 
current)  to  account  fully  for  the  enhancement  measured.® 

In  this  paper  we  describe  in  detail  a  new  and  quantitative  explanation^’1 
based  on  the  observation  that  current  enhancement  at  high  gas  pressures  is 
associated  with  large  transverse  excursions  of  the  beam.3-®  The  resistive 
hose  instability  which  drives  these  excursions  is  particularly  troublesome  at 
gas  pressures  above  a  few  torr,  and  develops  from  magnetic  interactions 
between  the  beam  current  and  plasma  eddy  currents  that  are  inductively 
generated  by  transverse  beam  motion.  Hosing  of  the  beam  spatially 
redistributes  the  net  current  density  and  thereby  alters  the  effective  system 
inductance.  A  rapid  drop  in  the  inductance  causes  the  monopole  electric  field 
(in  an  azimuthal  Fourier  expansion)  to  reverse  direction  and  drive  plasma 
current  parallel  to  the  beam  current.  Dipole  and  higher-order  fields 
complicate  the  redistribution  of  the  current  density  but  contribute  little  to 
the  spatially  integrated  net  current. 

According  to  the  present  model,  current  enhancement  is  thus  a  natural 
outgrowth  of  large-amplitude  hose  instability.  The  inductive  effects  which 
produce  the  instability  also  produce  current  enhancement.  The  enhancement  is, 
however,  a  nonlinear  effect  which  is  significant  only  for  large  beam 
displacements  and  cannot  be  seen  with  linearized  codes  and  models1 1-1 4  that 
are  valid  only  during  the  onset  and  initial  growth  of  the  hose  instability. 


\  % 


In  this  paper  we  present  analytic  and  numerical  treatments  of  inductive 
current  enhancement.  The  circuit  model  in  Sec.  2  provides  a  simple  yet  sound 
physical  basis  for  current  enhancement  due  either  to  beam  expansion  or  beam 
displacement.  In  Sec.  3  we  give  a  complete  analytic  field  solution  for  beam 
displacement  in  a  uniform  conductivity  channel.  Here  we  show  that  current 
gains  of  three  or  more  are  possible  if  the  beam  actually  strikes  a  metallic 
wall.  Consideration  of  beam  dynamics  and  conductivity  evolution  suggests, 
however,  that  the  hose  instability  is  unlikely  to  produce  current  gains  much 
above  a  factor  of  two.  In  Sec.  4  we  present  self-consistent  numerical 
simulations  of  beam  propagation  using  a  new  nonlinear  particle  code,  SARLAC, 
which  treats  both  large- amplitude  hose  instability  and  current  enhancement. 
Substantial  enhancement  (-  70  percent)  is  seen  for  beam  displacements  out  to 
half  the  wall  radiu3,  in  agreement  with  the  ETA  experiments.  Because  cause- 
and-effect  relationships  between  current  enhancement  and  various  beam  and 
cavity  parameters  are  obscured  by  the  complexity  and  sensitivity  of  the  hose 
instability,  we  also  present  simplified  simulations  using  a  prescribed 
dynamical  beam  displacement  to  further  elucidate  and  verify  the  analytic 
models. 


2.  Circuit  Analysis 


We  seek  in  this  section  to  provide  a  simple  physical  basis  for  current 
enhancement  due  either  to  beam  displacement  or  expansion.  A  few  models 
suffice  to  demonstrate  the  effect. 

A  circuit  model  is  often  used  to  compute  the  net  current  In  as  a  function 
of  beam  current  I&.  For  stable  beams  of  constant  radius,  the  circuit  model 
reduces  at  high  gas  pressure  to1 
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where  L  is  the  circuit  inductance,  R  is  the  plasma  resistance,  and  (In  -  Ib) 
is  the  plasma  current.  This  equation  relates  changes  in  the  magnetic  flux  LIn 
to  diffusion  of  magnetic  field  through  the  plasma  resistance  R.  An  important 
consequence  of  this  equation  is  that  the  maximum  In  is  less  than  the  maximum 
I.,  (provided  IR  <  I Initially).  Equation  (1)  therefore  precludes  current 
enhancement : 

F  .  !<«*)/!<■»)  < 
n  b  - 

Equation  (1  )  is  valid,  however,  only  if  the  inductance  L  is  constant.  In 

general,  the  inductance  is  not  constant  but  is  a  geometrical  parameter  which 

depends  on  the  spatial  distribution  of  the  beam,  plasma,  and  boundary 

currents.  The  boundary  current  is  presumed  to  reside  in  a  metallic  pipe  which 

encompasses  the  beam  and  plasma.  3eam  displacement  and  expansion  alter  the 

current  distribution  and  thereby  alter  the  inductance  L  and  flux  Ltn.  To 

3L 

include  these  effects  an  additional  term  proportional  to  I  -rr  must  be  added 

n  o  t 

to  the  left-hand  side  of  Sq.  (i).  This  additional  term  is  what  produces 
current  enhancement.  Such  a  possibility  has  been  previously  recognized  but 
5,15 


not  quantified. 


We  consider  first  the  case  where  the  plasma  conductivity  is  sufficiently 


high  to  suppress  electrostatic  effects  everywhere  inside  a  perfectly 
conducting  (metallic)  pipe  at  radius  b.  In  addition  we  make  the  long- 
wavelength  and  low-frequency  approximations: 
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where  c  is  the  speed  of  light.  Maxwell's  equations  then  reduce  to  the 
magnetostatic  Faraday-Ampere  law  for  the  axial  vector  potential  A,: 
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with  the  axial  electric  field  given  by 
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The  plasma  conductivity  <j  and  beam  current  density  are  specified  functions 
of  (r  ,  t).  The  boundary  condition  is 
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at  r  »  b.  In  subsection  2.3  we  extend  the  analysis  to  include  electrostatic 
effects  which  arise  if  the  medium  is  poorly  conducting  or  non-conducting  at 
radii  that  are  large  but  still  within  the  pipe.  A  complete  numerical 
treatment  is  given  in  Sec.  4. 


Equations  (3)-(5)  may  be  reduced  to  a  simple  circuit  equation  provided 
the  plasma  conductivity  is  azimuthally  symmetric  about  the  pipe  center.  To 


show  this,  consider  the  azimuthally  averaged  (monopole)  component  of  the 
electric  field: 
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3y  suitably  integrating  Eq.  (3)  and  applying  Sq.  (4)  and  boundary  condition 
(5),  one  finds  that 
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where  I  is  the  net  current  flowing  within  radius  r: 


rr  r2ir 
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The  plasma  current  equals  the  difference  between  the  net  current  and  the 
beam  current : 


ln  ~  rb  *  /  dr  r  J  d0  (9) 
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For  azimuthally  symmetric  conductivity,  only  the  monopole  field  contributes. 
Hence, 
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This  equation  or  its  equivalent  applies  to  ail  the  analytic  work  presented 
herein  but  not  to  the  numerical  simulations  in  Sec.  4  which  allow  the 


conductivity  to  be  azimuthally  asymmetric. 

To  express  Eq.  (10)  in  a  fora  similar  to  Eq.  (1),  we  define  the  radial 
profiles  of  the  net  current  and  plasma  conductance,  respectively,  by 


I(r , t)  -  I(r,t)/In(t) 


(11 ) 


and 
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where 


R(t)  «  ({  dr  2irraf1  (13) 
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is  the  plasma  resistance  per  unit  length.  Equation  (10)  may  then  be  rewritten 
as 
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where  we  define  L(t)  and  L(t)  as 
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and 
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Integration  by  parts  and  use  of  definition  ('.2)  produces  an  effective  system 
inductance  given  by 


LCt) 


and  a  parameter  L  given  by 


L(i) 
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Circuit  equation  (14)  together  with  definitions  (13).  (15),  and  06)  is 
an  equivalent  and  exact  representation  of  the  Faraday-Ampere  laws  (3)-(5), 
provided  that  the  plasma  conductivity  is  azimuthally  symmetric  about  the  pipe 
center.  Even  for  asymmetric  conductivity,  Eq.  (14)  is  often  a  good 
approximation.  In  this  equation,  the  effective  system  inductance  L  is  a 
geometrical  parameter  determined  by  the  location  of  the  pipe  boundary  and  by 
the  spatial  distributions  of  the  plasma  conductivity  and  net  current.  The 
parameter  L.,  which  was  absent  from  the  original  circuit  equation  (1), 


represents  a  change  in  inductance  due  to  redistr ibution  of  the  net  current. 


3L> 


Note  however  that  L  does  not  equal  unless  the  distribution  of  plasma 


A 

conductivity  is  constant  (-^  2,  *  0).  Current  enhancement  is  possible  only 


if  L  <  0. 


The  circuit  model  offers  a  simple  yet  potentially  accurate  means  for 
computing  the  net  current.  The  accuracy  is  limited  primarily  by  our  ability 
to  determine  R,  L,  and  L.  which  depend  (weakly)  on  the  spatial  distribution  of 
the  plasma  current  and  conductivity.  In  this  paper  we  make  only  partial  use 
of  the  circuit  equation.  We  nonetheless  provide  in  Appendix  A  some  ir.ductanc 
formulae  which  demonstrate  the  dependence  of  L  on  geometry.  When  using  these 
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formulae  to  compute  the  net  current  for  a  prescribed  beam  perturbation,  the 


parameter  L  should  generally  be  assigned  a  value  from  one-half  to  one 


times  §£. 


2.1  Enhancement  Due  to  Beam  Expansion 

To  demonstrate  the  utility  of  the  circuit  model,  consider  current 
enhancement  arising  from  self-similar  expansion  of  an  on-axis  beam  and 
plasma.  This  example  is  chosen  for  illustrative  purposes  only.  We  later 
conclude,  from  consideration  of  beam  dynamics,  that  beam  expansion  due  to 
axisymmetric  instabilities  such  as  sausage  and  hollowing  is  unlikely  to 
produce  current  gains  above  unity.  The  high  current  gains  observed 
experimentally  are  shown  in  subsection  2.2  and  following  to  arise  from  hose- 
induced  beam  displacement  with  expansion  as  a  secondary  effect. 

For  self-similar  expansion,  the  distribution  functions  for  current  and 
conductance  become  equal: 
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Equations  (15)  and  (16)  then  yield 
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so  that  circuit  equation  (14)  reduces  to 
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An  alternative  fora  is  the  energy  conservation  relationship 
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where  x  L  I  is  the  stored  aagnetic  field  energy  and  [(I  -  L  )R]  I  is  the 

2  n  no  n 

rate  at  which  the  current  I„  extracts  energy  ohaically  from  the  inductive 
electric  field  (In  -  Ib)R. 

Rapid  expansion  maximizes  current  enhancement  and  reduces  Eq.  (20)  to 
conservation  of  magnetic  field  energy: 

6 (L  In2)  -  0.  (21  ) 

Inserting  inductance  formula  (A1)  of  Appendix  A  into  this  relationship  yields 
a  maximum  current  gain  for  an  expanding  on-axis  beam  and  plasma  given  by 
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Here  aQ  i3  the  initial  beam  radius,  a  is  the  expanded  beam  radiu3, 

?  *  I  /I.  <  1  is  the  current  gain  immediately  or:or  to  expansion,  and 

o  no  3  “ 

2  * 

a  <  ^  is  a  geometrical  constant.  Examination  of  Eq.  (19)  reveals,  however, 
that  result  (22)  is  valid  only  if  the  beam  and  plasma  expand  in  a  time  short 


compared  with  the  monopole  decay  time, 


Otherwise,  ohmic  heating  extracts  sufficient  energy  during  the  expansion  that 
relationship  (21 )  is  a  poor  approximation. 


For  the  ETA  experiments  where  P/a  *  10,  Eq.  (22)  predicts  that  current 

o 

enhancements  of  30  percent  (F  »  1.3)  are  possible  if  F  *  1  and  if  the  beam 
radius  doubles  (a/a0  >  2)  in  less  than  a  monopole  decay  time  t  .  The 
simulation  results  presented  in  Sec.  4.1  show  such  expansion. 

Several  observations  should  be  made  regarding  current  enhancement  due  to 
beam  expansion.  The  first  is  that  axisymmetric  instabilities  such  as  the 
sausage  and  hollowing  modes1®-1®  can  cause  the  beam  to  expand  and  inductively 
raise  the  net  current.  Yet  such  expansion  is  unlikely  to  produce  current 
gains  above  unity.  The  reason  is  that  beam  expansion  due  to  the  axisymmetric 
resistive  instabilities  occurs  because  on-axis  plasma  current  initially  flows 
opposite  to  the  beam  current  and  magnetically  repels  it;  as  a  result,  Fq  <  1/2 
at  the  onset  of  the  axisymmetric  instabilities.1®'1®  Although  beam  expansion 
increases  F,  this  in  turn  strengthens  the  magnetic  pinch  force  on  the  beam 
which  restrains  or  reverses  further  expansion.  Usually* u  a  dynamic  balance  is 
reached,  or  further  radial  oscillations  occur,  but  with  F  <  1 . 

A  second  observation,  direct  from  result  (22),  is  that  the  pipe  radius  b 
is  the  radial  scale  length  governing  current  enhancement.  This  is  because  the 
inductance  L  changes  appreciably  only  when  the  beam  radius  (or  beam 
displacement,  as  shown  in  Sec.  2.2)  is  non-r.egligibie  compared  with  b.  One 
might  surmise  that  current  enhancement  cannot  occur  if  the  pipe  is  absent  or 
very  large.  This  conclusion  is  incorrect,  however,  as  shown  in  Sec.  2.3  where 
we  incorporate  electrostatic  effects  into  the  circuit  analysis.  There  we  show 
that  b  is  really  the  monopole  boundary  radius  at  which  the  as  i. mutually 
averaged  field  Ez0  falls  to  zero,  and  that  b  is  finite  if  the  medium  becomes 
non-conducting  at  large  radii  (ra  ♦  3  as  r  ■*  *). 


A  third  point  is  that  beam  expansion  is  a  secondary  effect  of  the 
resistive  hose  Instability:  as  the  beam  displaces  far  off  axis,  it  separates 
from  the  centroid  of  the  net  current  so  that  the  pinch  force  on  the  beam 
weakens  and  the  beam  expands.  The  circumstances  of  this  expansion  are 
sufficiently  complex  that  Eq.  (22)  should  be  treated  as  but  a  rough  guide  for 
current  enhancement  due  to  beam  expansion,  as  distinct  from  enhancement  due  to 
beam  displacement  alone. 

A  fourth  observation  is  that  the  energy  source  responsible  for  producing 
and  driving  the  plasma  currents  is  not  kinetic  beam  energy  directly  but 
magnetic  field  energy.  In  the  case  of  beam  expansion  the  stored  magnetic 
field  energy  not  only  drives  the  plasma  currents  but  also  accelerates  the  beam 
to  higher  energy. 

2.2  Enhancement  Due  to  Beam  Displacement 

Consider  now  current  enhancement  resulting  from  transverse  displacement 
of  a  rigid-rod  beam  of  fixed  radius  a.  We  again  assume  that  the  plasma 
conductivity  is  azimuthally  symmetric  about  the  pipe  center.  Only  the 
monopole  electric  field  then  contributes  to  the  net  current.  For  other 
conductivity  configurations  the  dipole  and  higher-order  fields  contribute  as 
well . 

A  simple  model  which  provides  a  semi-quantitative  estimate  of  current 
enhancement  i3  as  follows.  Imagine  a  beam  which  originally  resides  in  a  fixed 
plasma  channel  centered  within  a  pipe  of  radius  b.  The  plasma  channel  and 
beam  have  a  small  radius  a  <<  b. 

A  sudden  displacement  of  the  beam  to  an  off-axis  location  r&  perturbs  the 
plasma  currents  and  wall  currents.  The  values  of  these  currents  immediately 
after  the  displacement  can  be  readily  computed  provided:  (i)  the  displacement 


occurs  instantaneously;  and  (ii)  the  displaced  beam  is  wholly  outside  the 
plasma  but  wholly  inside  the  pipe,  i.e., 


2a  <  r.<  b  -  a. 
b 


The  plasma  then  behaves  as  a-  perfect  conductor  while  the  beam  can  be  treated 

as  a  line  current  at  r.  . 

-b 

A  circuit  treatment  of  this  problem  is  given  in  Appendix  A.  Here  we  give 
for  comparison  a  field  treatment.  We  begin  by  dividing  the  plasma  current 
into  two  components.  The  homogeneous  component  is  derived  by  imagining  that 
the  beam  current  1^,  is  suddenly  moved  from  inside  the  plasma  to  outside  the 
pipe.  Because  magnetic  flux  remains  trapped  within  the  plasma,  the 
homogeneous  plasma  current  is  given  by 


Iph  '  ‘no 


where  Iro  is  the  net  current  Just  prior  to  beam  displacement. 

The  driven  component  of  the  plasma  current  is  derived  by  imagining  that 
the  line  current  Ip  is  suddenly  moved  from  outside  the  pipe  to  the  location 
inside.  For  thi3  problem  the  vector  potential  Az  is  zero  both  at  the  pipe 
boundary,  r  *  b,  and  everywhere  within  the  plasma,  r  £  a.  These  two  boundary 
conditions  determine  the  driven  (surface)  plasma  current  I  .. 

pu 

We  evaluate  Ipd  by  integrating  the  Faraday-Ampere  law  (3)  over  r  and  3  to 
obtain  for  r  >  a: 


dF2  •  -  -f  H<r  •  V! 

C  P 
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where  H  is  the  Heaviside  step  function  and  where  the  monopole  vector  potential 
is  defined  by 


2ir 


Vr)  ■  k  /  de  Az(r'9) 

o 


(27) 


Integrating  Eq.  (26)  and  applying  the  boundary  conditior 


A  (a)  -  0 
zo 


(28) 


yields 


A  (r)  ■  “  ln(f*/a)  ♦  L  H(r  -  r  )  ln(r/r  )]. 

zo  2  pd  b  b  b 

c 


(29) 


Applying  the  boundary  condition 


A  (b)  -  0 
zo 


(30) 


produces 


I.  ln( b/r  ) 

pd  _ b_ 

I.  ln(b/a)  * 

b 


(31 ) 


Solution  (31 )  for  the  driven  plasma  current  I ^  has  a  simple  physical 
basis.  Injecting  a  current  source  between  two  conducting  surfaces  drives 
return  current  in  each  of  the  surfaces.  The  total  return  current  equals 
-Ib.  The  fraction  of  return  current  flowing  in  either  one  of  the  surfaces  is 
a  function  of  the  surface  geometry  and  of  the  source  location.  For  coaxial 
conductors  of  inner  ar.d  outer  radii  a  and  b,  the  fraction  flowing  on  the  inner 
(plasma)  conductor  is  giver,  by  Eq.  (31). 
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'  l  l  -l'JV  P  V  * 


The  overall  net  current  inside  radius  b  is  the  sun  of  the  beam  current 
and  the  plasma  currents  Iph  and  Ip<1.  From  Sqs.  (25)  and  (31)  we  thus  have 


I  +  I  *  I 
b  ph  pd 


F_  ♦ 


lr.(rb/a) 


0  lr.(b/a) 


(32) 


where  ?3  -  Ino/Ib  is  the  current  gain  immediately  prior  to  pushing  the  source 

off  axis.  Typically  ?0  <  1  due  to  plasma  return  currents  which  were  induced 

when  the  or.-axis  source  was  originally  injected  into  the  plasma. 

Result  (32)  for  the  current  gain  F  was  derived  assuming  instantaneous 

beam  displacement.  We  show  in  Appendix  A,  however,  that  magnetic  flux  relaxes 

slowly  on  a  monopole  decay  time,  t  -  L/R.  Result  (32)  can  thus  3till  be 

0 

applied  provided  the  beam  moves  off  axis  to  r&  in  a  finite  but  short  time 

i(rj  <  [ln(b/a)  *  a].  (33) 

~b  0  Rc2 

Although  the  net  current  and  monopole  fields  relax  slowly,  the  dipole  and 
higher-order  fields  (in  an  azimuthal  Fourier  decomposition)  relax  rapidly  on  a 
dipole  decay  time  given  by1 1 

-  (2RC2)"1 .  (3M 

Comparison  of  Eqs.  (33)  and  (34)  reveals  that  the  mor.opole  and  dipole  decay 
times  are  typically  related  by 

VTd  "  I*Cln(b/a)  ♦  a]  »  1.  (35) 


The  plasma  current  thus  spatially  redistributes  before  it  decays. 
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The  distinction  between  decay  and  redistribution  of  the  plasma  eddy 
currents  is  important  in  understanding  hose-induced  current  enhancement.  The 
resistive  hose  instability  develops  from  a  phase  lag  between  dynamical  beam 
displacement  and  redistribution  of  the  dipole  fields.  In  the  absence  of  a 
phase  lag,  the  dipole  fields  produce  a  restoring  force  which  causes 
oscillation  but  no  growth  in  beam  displacement. 1 ^  The  dipole  decay  time 
rd  controls  the  phase  lag  and  thus  characterizes  the  growth  rate  of  the 
instability.  During  a  monopole  time  tQ  >>  xd,  large-amplitude  hose 
displacement  is  possible.  A  hose-unstable  beam  of  long  duration  <»  V  13 
thus  likely  to  produce,  given  sufficient  propagation  distance,  a  significant 
increase  in  the  net  current. 

For  the  ETA  experiments  where  b/a  *  10,  Eq.  (32)  predicts  that  current 
enhancements  of  30  percent  (F  -  1.3)  are  possible  even  for  modest 
displacements  of  a  few  beam  radii,  rfa  •  2a.  A  current  enhancement  of  nearly 
100  percent  (F  +  2)  is  possible  for  large  displacements,  r&  ♦  b  -  a. 

2.3  Seams  Propagating  in  Unbounded  Neutral  Gas 

In  the  preceding  treatments  we  assumed  that  a  metallic  drift  tube  was 
present  and  that  the  gas  conductivity  everywhere  within  the  drift  tube  was 
sufficiently  large  to  ensure  space-charge  neutralization.  If  this  is  not  the 
case,  electrostatic  terms  must  be  included  in  the  field  and  circuit  equations. 
We  show  here,  using  a  revised  circuit  analysis,  that  the  phenomenon  of  hose- 
induced  current  enhancement  is  qualitatively  unchanged  even  if  the  drift  tube 
is  absent  or  lies  far  outside  the  region  where  space-charge  neutralization 
occurs.  The  principal  revision  is  that  the  monooole  boundary  radius  b,  where 
Ezo  ■  0,  is  not  the  pipe  radius  but  a  "vacuum  radius"  beyond  which  the  space- 
charge  fields  are  unsuppressed. 


«.  vie* 
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The  revised  circuit  analysis  is  based  on  Maxwell's  equation  for  the  axial 


electric  field: 


V2  E 
1  2 


,2- 

d  e. 


3z‘ 


c2  at2 


It  (Jb  +  aEz)  +  4ir  If 

c 


(36) 


55 

where  p  is  the  charge  density.  As  shown  by  E.  P.  Lee,  this  equation  can  be 
simplified  for  beams  which  are  both  ultrarelativistic, 


8  *  v  /c  -►  1 , 
2* 


(37) 


and  paraxial, 


a  d 
Be  dt 


0(a/A  >  «  1. 


(38) 


Here  v2  is  the  beam  axial  velocity,  a  is  a  characteristic  beam  radius, 
t  *  t  -  z/3c  is  a  temporal  variable  designating  a  particular  beam  slice, 

■a  i  /2 

and  X  *  2ira(8Ymc  /el  )  is  a  characteristic  betatron  wavelength.  The 

3  »• 

2  -1/2 

electron  charge-to-mass  ratio  is  e/m,  and  Y  *  (1  -3  )  is  the  usual 
relativistic  factor. 

3y  changing  variables  from  (z,t)  to  (z,t)  and  by  setting  3  -  1  and 
(t-J  -  0,  we  can  rewrite  Eq.  (36)  as 

oZ  T 


72e  •  |-  (J 

1  z  2  3t  b 
c 


oE  -  pc) 
z 


(39) 


-? 

which  is  typically  valid  to  order  Y  <<  1  and  a/X.  <<  1 .  The  neglect  of 
axial  derivatives  is  denoted  the  "frozen-field  approximation"  and  follows  from 
the  paraxial  condition  (38).  The  full  Lee  field  equations,  which  are 
presented  in  Sec.  4,  employ  an  additional  approximation  not  needed  here. 


For  azimuthally  symmetric  conductivity,  we  may  duplicate  the  procedure  in 


Sec.  2  and  reduce  Eq.  (39)  to  a  circuit  equation  for  the  net  conduction 
current  I„: 


L 


3_ 

3t 


I 

r. 


♦  LI 


(In  -  VR 


(40) 


where  the  effective  system  inductance  is  now  defined  by 


ut>  -  %/f  i  <i-i) 

c  o 


(41) 


and  the  parameter  L  is  defined  by 


C<T)  ■  h  S  f  I  It  (I  - 


c  o 


Here  the  radial  profile  of  the  distributed  charge  is  given  by 


2ir 


Qfr.r)  -  9  1  |  dr  r  J  d0 


(42) 


(43) 


while  the  distributed  net  charge  is  given  from  charge  conservation  by 


2ir 


Qn(x)  -  /  dr  r  J  de  p  »  In(T)/c. 


(44) 


The  revised  circuit  equations  (uo)-(42)  closely  resemble  the  earlier  versions 
( i 4 )-( i 6 ) .  The  difference  is  the  presence  of  Q  and  the  apparent  absence  of  a 
finite  cutoff  radius  b  in  the  inductance  formulae.  Recall  that  a  finite 


cutoff  is  required  for  current  enhancement  to  occur. 


A  finite  cutoff  in  fact  exists.  At  some  large  radius  b  out3ide  the  beam, 
the  conductivity  becomes  so  low  that  the  medium  act3  a3  a  vacuum:  typically, 

4ir<jr/c  <0.1  (45 ) 

A  A 

for  all  r  ^  b.  3eyor.d  the  vacuum  radius  b,  the  profile  functions  I  and  Q 
equal  unity  by  definition.  Inductance  formulae  (41 )  and  (42)  can  thus  be 
terminated  without  approximation  at  r  -  b. 

For  intense  tightly  pinched  beams,  conductivity  in  the  corona  region 
outside  the  beam  ia  generated  both  from  collisional  ionization  by  the  expanded 
beam  head  and  from  avalanche  ionization  in  the  radial  electrostatic  field. 

For  such  beams,  the  vacuum  radius  b  is  typically  large  (>>  a)  and  insensitive 
to  displacement  or  expansion  of  the  beam  body.  Moreover,  the  corona 
conductivity  rises  rapidly  inside  b  so  that  space  charge  accumulates  near 

A 

r  -  b.  Behind  the  beam  head,  the  charge  density  p  and  profile  Q  therefore 
approach  zero  for  r  <  b,  and  Eqs.  (40)-(42)  become  virtually  indistinguishable 
from  Eqs.  (I11)- (16).  Current  enhancement  for  a  given  perturbation  of  the  beam 
body,  where  large  hose  growth  is  most  probable,  is  thus  almost  the  same 
whether  or  not  a  metallic  boundary  is  present  at  r  -  b  (or  at  r  >  b). 

The  presence  or  absence  of  a  metallic  wall  beyond  the  vacuum  radius  does 
not  directly  affect  current  enhancement  because  a  vacuum  boundary,  like  a 
metallic  boundary,  shorts  out  the  monopole  electric  field  E,0  which  determines 
the  net  current  for  axisymmetric  conductivity.  Sven  for  asymmetric 
conductivity,  the  monopole  contribution  usually  dominates  because:  (i)  o  is 
positive  definite  and  thus  has  a  large  monopole  component;  and  (ii)  the  dipole 
and  higher-order  azimuthal  components  of  E2  decay  much  more  rapidly  than 
Ezo.  Circuit  equations  (4Q)-(42)  are  thus  usually  a  fair  approximation  to  Eq. 
(39)  even  if  <j  is  asymmetric. 
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The  nature  of  the  boundary  can,  of  course,  strongly  affect  the  dipole  and 
higher-order  fields  which  control  beam  displacement.  Large  displacements  are 
restrained,  for  example,  by  metallic  boundaries  which  carry  a  return  current 
that  magnetically  repels  the  beam,  but  are  not  restrained  by  vacuum  boundaries 
where  the  return  current  and  associated  wall  forces  are  absent  or  weak.  The 
nature  of  the  boundary  can  thus  alter  the  degree  of  hose  growth  and  consequent 
current  enhancement,  but  does  not  alter  the  occurrence  or  qualitative  behavior 


of  the  phenomenon. 


3.  9eam  Displacement  in  a  Uniform  Plasma 

In  this  section  we  consider  a  broad  channel  with  uniform  conductivity 
extending  out  to  the  metallic  boundary.  The  extension  in  a  has  a  pronounced 
effect  on  current  enhancement  and  represents  the  opposite  limit  to  the  narrow 
conductivity  channels  considered  in  Sec.  2.2.  Unfortunately,  a  circuit 
analysis  is  impractical  for  extended  o  because  the  spatial  distribution  of  the 
plasma  current  can  no  longer  be  estimated  a  priori.  We  give  instead  a 
complete  field  solution  for  a  problem  which  is  analytically  tractable:  viz., 
an  infinitesimally  thin  current  source  moving  at  constant  transverse  velocity 
u  through  a  medium  of  uniform  and  constant  conductivity  o.  The  underlined 
words  represent  the  simplifying  assumptions.  In  subsection  3.2  we  discuss  the 
consequences  of  relaxing  some  of  these  assumptions. 

The  solution  to  our  model  problem  can  be  expressed  in  terms  of  two 

dimensionless  parameters:  r&(t)/b  and  5Q/b  where  rb(t)  is  the  beam 

2 

displacement  and  <SQ  *  c  /4ircu  is  identified  as  a  skin  depth.  The  solution 
would  depend  on  a  third  parameter,  a/b,  if  the  source  were  of  finite 
thickness. 

Solutions  are  given  with  and  without  displacement  currents  included  in 
the  analysis.  As  expected,  the  displacement  currents  are  unimportant  provided 
u/c  <  0.1.  We  show  in  Appendix  C  that  the  solutions  are  insensitive  to  the 
cavity  shape  as  well. 

We  describe  this  problem  and  its  solution  in  considerable  detail  because 
it  provides  additional  insight  and  because  it  reveals  two  new  features.  The 
first  is  that  current  enhancement  in  broad  conductivity  channels  does  not 
become  appreciable  until  the  current  source  approaches  to  within  a  skin  depth 
of  the  conducting  boundary.  The  3ecor.d  and  more  surprising  feature  is  that 
for  o  /b  <<  1 ,  a  current  gain  of  three  or  more  is  attained  as  the  source  is 


about  to  strike  the  boundary.  Recall  that  Eq.  (32)  restricts  the  current  gain 
for  a  narrow  (and  stationary)  conductivity  channel  to  a  maximum  value  of  two. 
Extending  the  plasma  conductivity  out  to  the  boundary  thus  reduces  the  current 
gain  for  modest  displacements,  r  <  b/2,  but  increases  the  gain  for  large 

D 

displacements  r^  -*■  b. 

3.1  Formulation  and  Solution 

For  simplicity  we  consider  a  one-dimensional  system  in  slab  geometry. 

Ilhe  extension  to  a  two-dimensional  system  may  be  found  in  Appendix  C.]  In 
our  model,  an  infinitesimally  thin  current  sheet  of  surface  density  Ik  flows 
in  the  z-direction.  The  current  sheet  is  initially  located  midway  between  two 
perfectly  conducting  plates,  ?L  and  PR,  at  x  -  0  and  x  »  1,  respectively.  We 
assume  that  the  medium  between  P^  and  Pg  is  characterized  by  a  constant 
electrical  conductivity  o.  We  further  assume  that  the  current  3heet  has  beer, 
situated  at  x  *  «  1/2  for  a  long  time  so  that  all  plasma  current  has 

decayed  to  zero  at  t  *  0.  For  t  Z  0,  the  current  sheet  is  given  a  uniform 
x-ward  motion,  i.e.,  x*j(t)  *  1/2  *  ut  where  the  speed  u  is  assumed  to  be 
constant. 

The  fields  are  governed  by  the  diffusion  equation  for  the  longitudinal 


The  initial  condition  is 


A  (x,0)  »  (irilb/c)  (1  -  1 2x  -  l\/l). 


(48) 


reflecting  the  assumption  that  a  steady  magnetic  field  is  established  prior  to 
t  -  0. 

The  current  enhancement  factor  is  given  by 

F(t)  -  1  -  («/I.e)  J  dx  3A  /3t,  (49a) 

b  o  2 


where  -(a/c)3A  /3t  is  the  current  density  induced  within  the  conductivity 
z 

channel.  Alternatively, 


F(t) 


,  3A 

]_  r  e _ z 

lw  L  4ir  3x 


x*i 


3A 
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4ir  3x  1  J 


1  x=0 


(49b) 


where  the  square  bracket  represents  the  wall  currents.  Equation  (49b)  is 
readily  obtained  by  integrating  Eq.  (46)  from  x  *  0  to  x  *  1. 

The  diffusion  equation  (46),  together  with  its  associated  initial 
conditions  and  boundary  conditions,  is  solved  in  Appendix  3.  There  we  show 
that  just  before  the  current  3heet  makes  contact  with  P^,  the  current 
enhancement  factor  is  given  by  [cf.  Eq.  (313)] 


:ar.h 


rTV-i 
'2  j 


J3 


(50) 


where 


v 


4uu£/c^ 


(51  ) 
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is  the  ratio  of  the  characteristic  diffusion  time  x 


*  2oZ  /c  to  the 


dif 

transit  time  x,.  *  l/2u.  The  infinite  sum  in  Eq.  (50)  converges  rapidly  for 
finite  values  of  v. 

The  enhancement  factor  F,  as  given  by  Eq.  (50),  is  plotted  in  Fig.  1  as  a 
function  of  the  normalized  velocity  v.  The  most  interesting  and  unexpected 
feature  is  that  F  approaches  an  asymptotic  limit  of  three  as  v  becomes  large. 
From  Fig.  1,  one  sees  that  F  >  2  for  v  ^  1.  For  v  <  0.5,  a  good  approximation 
is  F  »  1  +  nv/2,  as  is  readily  verified  from  Eq.  (50)  and  is  also  easily 
deduced  from  a  quasi-static  argument. 

The  time  evolution  of  F  gives  additional  insight  into  the  physical 
processes.  We  show  in  Appendix  B  that 


1  ♦  JI 

IT  . 

n*1 
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1-1)  /  2  rnitt'i  .  rnirti 

- 5 - j-  {-v  cos(— J  x  nv  sin(— j 

n(n  x  v2)  2lt  2Tt 


2 

+  v  e 
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This  infinite  series  is  absolutely  and  uniformly  convergent  for  all  t  ^  0. 

The  evolution  of  F  as  a  function  of  normalized  time  t/x,.  is  shown  in  Fig.  2 
for  various  values  of  v.  Observe  that  current  enhancement  becomes  appreciable 
only  when  x^  reaches  within  a  "boundary  layer"  width  of  PR,  i.e.,  when 
t  -  xk  <  x.  /v.  The  width  of  this  boundary  layer  (skin  depth)  is  on  the  order 
of  l/v,  as  suggested  by  Fig.  2.  This  may  also  be  deduced  by  a  dimensional 
argument  from  the  diffusion  Eq.  (d6). 

In  Appendix  C  we  extend  our  model  to  two  transverse  dimensions  and 
represent  the  beam  as  a  current-carrying  wire  which  moves  at  constant 
transverse  velocity  inside  a  waveguide  filled  with  conducting  medium.  We  find 
that  all  qualitative  features  of  the  one-dimensional  model,  including  the 


maximum  enhancement  factor  of  three,  are  essentially  unchanged.  Indeed,  these 


features  are  of  such  a  universal  nature  that  they  are  independent  of  the 
initial  location  of  the  current  source  as  long  as  tdiI,  >  >  t^.  In  the  next 
subsection  we  discuss  the  effects  of  finite  beam  thickness  and  of  nonuniform 
velocity. 


3.2  Interpretation  and  Generalizations 

To  elucidate  the  preceding  unexpected  features  of  a  boundary-layer  effect 
and  F  -*•  3  for  v  >>  1 ,  we  decompose  the  induced  plasma  current  into  two  parts. 
The  first  part  is  associated  with  relaxation  of  the  magnetic  field  which  was 
present  prior  to  t  -  0.  This  part  is  independent  of  the  motion  of  the  current 
source.  The  second  part  is  due  to  the  transverse  motion  of  the  current 
source,  but  is  independent  of  the  initial  magnetic  field.  Such  a 
decomposition  also  enables  us  to  predict  the  qualitative  behavior  when  the 
beam  motion  is  nonuniform  or  the  beam  size  is  finite.  These  generalizations, 
together  with  the  influence  of  the  displacement  currents,  will  be  addressed 
later. 

The  above  mentioned  decomposition  is  equivalent  to  expressing  the 
solution  Az(x,t)  to  Sq.  (46)  as  a  superposition: 


Az(x,t)  =■  Ah(x,t)  *  Ad(x,t). 


(53) 


In  Eq.  (53).  A^(x,t)  is  the  homogeneous  solution  and  Ad(x,t)  is  the  driven 
solution.  We  similarly  decompose  the  plasma  current  density  as 
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^  ^  1  all  <  ^  a. .!  m-f  1.1.  t-A.I  »  »  -  *  -  »  -  *  -  *  -  "  ^  ^  ^ 


and  the  plasma  current  as 
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(55) 


The  homogeneous  solution  A^  describes  the  relaxation  of  the  magnetic 
field  which  pre-exists  before  t  -  0.  It  is  governed  by 
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and  is  subject  to  the  boundary  condition 


Ah(°,t)  -  Ah(l,t)  -  0 

and  to  the  initial  condition 


(57) 


Ah(x,0)  -  (irZIb/c)  (1-  |2x-l|/l). 


(58) 


The  driven  solution  Ad  represents  the  response  to  the  transverse  motion 
of  a  current  source  which  is  "switched  on"  instantaneously  at  time  t  *  0. 
That  is,  the  initial  magnetic  field  plays  no  role  in  the  driven  solution  Ad. 
Thus  Ad  satisfies 
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and  is  subject  to  the  boundary  condition 


Ad(0,t)  -  AdU,t)  -  0 


(60) 


vv-Sv-i- 


26 


and  to  the  initial  condition 


Ad(x,0)  -  0. 


(61 ) 


Both  Ah  and  A^  may  be  obtained  by  the  method  outlined  in  Appendix  B.  Let 
us  first  concentrate  on  Ah.  It  is  given  by  [cf.  Eq.  (B15)3 


411. 
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2xJ  sinl-y-J 


(62) 


which  relaxes  on  the  time  constant  Ah  also  gives  rise  to  Jph  which 

decays  on  the  same  time  scale.  The  evolution  of  Jph  is  shown  in  Fig.  3  for 

two  values  of  v.  Its  spatial  distribution  is  approximately  gaussian.  Note 

that  for  t  <  <  tdift  I  •*  Ip,  as  expected  on  physical  grounds  and  readily 

verified  from  solution  (62).  Thus,  in  the  limit  v  >  >  1,  I  .  contributes  one 

ph 

unit  to  the  current  enhancement  factor  F. 

The  driven  solution  Ad  gives  rise  to  the  unexpected  features  mentioned  at 
the  beginning  of  this  subsection.  It  is  given  by 
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where 
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(64) 


with  h  *  it(1  *  t/T J/2  and  h  -  t/2.  The  evolution  of  J  .(x,t)  is  shown  i: 
to  pd 

Fig.  4  for  v  -  0.5  ar.d  for  v  -  4. 


We  note  that  J  .(x,t)  is  negative  for  t/r  <  <  1  but  becomes  positive 

PC  U 

for  t/rt  *  1 .  This  behavior  may  be  understood  by  observing  that  Jpd  is  the 
plasma  current  due  to  a  current  source  that  is  switched  on  at  time  t  »  0. 
Initially,  Ad  -  0.  An  induced  field  and  current  are  thus  generated  at  time 
t  «  0+  so  as  to  maintain  zero  magnetic  field  by  Lenz's  law.  Initially,  Ipd 
exactly  cancels  Ife,  i.e.,  I  d  -*•  - 1 p  as  t  *  0  .  A  major  portion  of  the 
negative  inductive  electric  field  lies  ahead  of  the  sheet  current,  especially 
when  v  >  >  1  [cf.  Fig.  4],  However,  as  this  inductive  electric  field  "wave 
front"  approaches  the  conducting  boundary  Ppj,  it  is  reflected  and  its  polarity 
is  changed  in  the  same  manner  as  when  a  voltage  wave  is  reflected  by  a  short 
on  the  path  of  a  transmission  line.  In  this  analogy,  the  transmission  line  i3 


lossy.  But  in  the  limit  of  high  electrical  conductivity  (v  >  >  1),  the 


reflected  current  pulse  has  not  decayed  substantially,  i.e.,  I  *  +  Ip  as 

t  •»  Tt  for  v  >  >  1 .  Hence,  F  -*  1  as  t  ■*  0,  while  F  ♦  3  as  t  ♦  xt  for  v  >  >  1 . 

In  fact,  a  simple  analytic  solution  for  the  driven  solution  Ad  and  the 
associated  current  density  Jpd  may  be  derived  in  the  limit  that  v  >  >  i . 

First,  it  is  easily  seen  that  Eq.  (59)  admits  a  wave-like  solution  of  the  form 
Ad(x,t)  <*  f  [x  -  xp(t)].  Such  a  solution  yields  an  induced  current  density 
given  by 


Jpd(x,t) 


x  <  xt 


•^b^i  -  irv(x-x.  )/Z 
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which  is  indeed  an  approximate  solution  when  the  current  sheet  is  far  from  the 

wall.  The  important  feature  is  that  Jp^  lies  entirely  ahead  of  the  beam  and 

thus  reaches  the  wall  before  the  beam.  Because  Ap(x,t)  does  not  satisfy  the 

boundary  condition  A^(Z,t)  »  0,  an  image  current  source  of  the  form 

-I.  5 { x-[2Z-x.  (t ) ] }  must  be  added  to  render  the  vector  ootential  zero  at 
0  0 


x  *  l.  Associated  with  this  image  is  the  plasma  current  density 


Jod(x*t) 


0  ,  x  >  21  -  x. 

D 

(1^)  e’"x  -  21  '  S)/l,  x  <  21  -  xs 


(66) 


which  represents  the  reflected  wave  previously  mentioned.  Summing  all 

1  2 

currents  (1^,  1^,  I  ^  Ip^)  between  P^  and  PR  yields  a  current  enhancement 

factor  of  3  as  xb  *  2,  since  the  contribution  from  I  ^  becomes  negligibly 

2 

small  while  the  contribution  from  I  becomes  significant.  Equations  (65)  and 
(66)  demonstrate  why  current  enhancement  occurs  within  a  boundary  layer  of 
width  2/v  of  PR. 

For  simplicity  of  exposition,  the  calculations  of  this  section  have  been 

based  on  idealized  model  assumptions.  We  now  relax  some  of  these  assumptions 

and  use  the  insight  provided  by  our  results  to  predict  what  happens  when  the 

beam  has  finite  thickness  or  nonuniform  velocity  to  the  wall.  (One  could  also 

draw  insight  for  the  case  of  non-uniform  conductivity.)  Consider  fir3t  the 

neglect  of  finite  beam  thickness.  Thick  beams  can  be  treated  as  a  linear 

superposition  of  thin  beams.  X  beam  is  thin  if  its  thickness  2a  is  small 

2 

compared  with  the  skin  depth  <5q  *  l/rtv  -  c  /Uirou.  If  the  beam  is  thick,  the 
limiting  value  of  F  is  reduced  because  part  of  the  beam  is  more  than  one  skin 
depth  away  from  the  wall  when  the  beam's  leading  edge  reaches  the  wall.  High 
current  enhancement  is  thus  possible  in  broad  conductivity  channels  only  if 


a  <<  c2/4itcu  <<  2. 


(67) 


Current  gains  above  three  are  possible  in  broad  channels  if  the  beam 
accelerates  to  the  wall.  To  show  this,  consider  a  case  where  xfc  maintains 
uniform  motion  at  speed  u  until  t  *  Tj  <  t^,  but  undergoes  an  impulse 
acceleration  so  that  dx./dt  -  u,  »  constant  (>  u)  for  t  >  T. .  For  time  t  < 

b  i  1 

T1t  the  solution  remains  the  same  as  that  obtained  in  the  previous  subsection. 


Specifically,  F(T.)  >  1.  For  time  t  >  T. , 


since  u.  is  uniform,  we  may  again 


use  the  previous  section;  the  only  modification  is  that  the  solution  is  time' 


iv 

'! 
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translated  by  T1  and  that  the  "initial  condition"  at  t  «  T,  must  reflect  a 
larger  initial  current,  F(T,)Ib.  Now,  for  u1  ♦  •,  we  have  as  xb  *  1 , 


F  -  (I.  ♦  I  .  I  ,  )/I.  -  1  ♦  1  ♦  (I_/IK)  •  2  *  F(T, )  >  3.  Here,  we  have 
b  pd  ph  b  ph  b  i 


used  the  result  that  I  .  -*  I.  irrespective  of  the  "initial  condition".  In 

pd  b 


retrospect,  it  is  somewhat  surprising  that  F  ■*  3  regardless  of  the  magnitude  of 
u  so  long  as  it  is  large  and  constant,  but  that  F  may  exceed  3  if  du/dt  >  0. 

We  saw  in  Sec.  2.2  that  the  current  gain  for  narrow  and  stationary 
conductivity  channels  is  restricted  to  F  <  2;  this  low  value  of  Fmax  occurs 
because  only  the  beam  current  1^  physically  displaces  to  lower  the  inductance 


L.  For  the  broad  channels  discussed  in  this  section,  Fmax  can  be  as  large  as 


three  or  more  because  a  portion  of  the  plasma  current,  Ipd,  moves  with  Ib  to 
further  lower  L.  Maximum  current  gain  is  achieved  if  all  current  moves  with 


so  that  the  final  value  of  L(>^  0)  is  minimized.  This  is  the  case  for  a 


narrow  conductivity  channel  which  physically  displaces  with  the  beam.  In  this 
case,  the  final  values  of  L  and  F  are  limited  only  by  the  finite  breadth  of 
the  beam  and  channel,  which  prevents  all  the  net  current  from  simultaneously 
reaching  the  boundary.  By  applying  flux  conservation  to  inductance  formula 
(A4)  of  Appendix  A,  we  conclude  that  in  cylindrical  geometry  the  maximum 
current  gain  is  given  by 


Zn(b/a)  *  a 


max 


ln(2)  *  a 


(63a) 


where  a  is  a  characteristic  beam-channel  radius,  b  is  the  boundary  radius, 


2  1 

a  <  — ,  and  the  beam  displacement  rp  varies  from  zero  to  vb  -  a).  Similar 


arguments  applied  in  planar  geometry  produce 


F  *  1/ a 
max 


(58b) 
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where  l  is  the  plate  separation  and  2a  characterizes  the  thickness  of  the  beam 
and  conductivity  sheets.  Hose-unstable  beams  are  unlikely,  however,  to 
produce  current  gains  much  above  a  factor  of  two  because  they  leave  behind  a 
highly  conducting  on-axis  channel  that  traps  plasma  current  and  magnetic  flux; 
see  Sec.  4. t . 

The  displacement  current  has  been  neglected  thus  far,  and  in  the 

numerical  simulations  of  the  next  section  as  well.  While  the  physical 

description  given  here  indicates  that  the  displacement  current  effects  should 

be  small,  we  nevertheless  repeated  our  calculations  with  the  displacement 
-12  2 

current  term  (Uirc)  3  A  /3t  included  on  the  left-hand  side  of  Eq.  (46).  In 

z 

addition  to  conditions  (47)  and  (48),  the  initial  condition  3A  /3t  -  0  at 

z 

t  ■  0  is  introduced  to  completely  specify  the  solution  Az.  The  current 
enhancement  factor  F  is  calculated  according  to  Eq.  (49b)  which  now  implicitly 
includes  the  displacement  current.  The  result  is  shown  in  Fig.  5.  We 
conclude  from  Fig.  5  that  the  displacement  current  generally  increases  the 
enhancement  factor  F,  but  that  F  does  not  exceed  a  value  of  3.2  if  u/c  <  0.3. 
This  explicit  demonstration  of  the  unimportance  of  the  displacement  current 
provides  additional  confidence  in  Lee's  simplified  field  equations  which  are 
used  in  the  simulations  presented  in  Sec.  4. 


J*.  Computational  Models 


The  resistive  hose  instability  has  been  widely  studied  using  linearized 

simulation  models1 14 * 1  20  which  are  valid  only  for  small  beam  displacements. 

However,  the  current  multiplication  effect  described  in  the  previous  sections 

is  an  inherently  nonlinear  phenomenon.  In  this  section,  we  describe  nonlinear 

simulations  using  a  new  particle  code,  SARLAC.  This  code  provides  a  self- 

consistent  treatment  of  the  hose  instability  for  large  displacements  and 

reproduces  the  current  enhancement  effect.  A  more  complete  description  of  tne 

21 

code  is  given  elsewhere. 

SARLAC  employs  a  number  of  approximations  which  have  been  widely  used  in 
studies  of  electron  beam  propagation.  Longitudinal  and  transverse  beam 
dynamics  are  decoupled  and  treated  separately  (to  order  v  /v  *  2ira/\  <<  1) 

Z  3 

by  invoking  the  ultrarelativist ic  approximation  (37)  and  the  paraxial 
condition  (38).  The  field  solver  is  based  on  the  Lee  field  equations:22 

72(a  *  *)  -  -  £  (J  -  o  f£)  (69) 

c  b  oc 


and 


72  M 

1  3? 


(70) 


which  are  obtained  by  dropping  derivatives  of  the  transverse  vector  potential 

f-r—  A  ■»  oj  and  by  employing  the  frozen-field  aoproximatior.  discussed  in  Sec. 

2.3-  Here  $  is  the  electrostatic  potential,  and  A  *  A  -  $.  The  independent 

z 

temporal  variable  used  in  SARLAC  is  not  x  but  its  equivalent,  z  *  cx  »  ct  -  z, 
which  is  the  distance  behind  the  beam  head.  The  axial  electric  field  is  given 


while  the  transverse  electric  and  magnetic  fields  are  given  respectively  by 


E  • 
'1 


(71b) 


and 


B  «  7  (A  +  $)  x  z. 

'1  ~± 

Metallic  boundary  conditions  are  assumed: 

A  -  *  -  0 


(7lo) 


(72) 


at  a  specified  pipe  radius  b. 

The  conductivity  a  is  treated  by  a  simple  rate  equation: 

-  KJ&  ♦  v^/c  -  0rpo2  (73) 

where  in  air  we  use  the  beam  collisionai  ionization  coefficient  K  -  5X1 0-1* 
cm/statamp-sec  and  the  recombination  coefficient  Sr  »  7X10  1 ^  sec/cm-atm.  The 
gas  density  p  is  specified  in  atm.  The  electron-avalanche  ionization  rate 
-  vf(p;E)  is  a  complicated  function  of  density  p  and  electric  field  E,  and 
is  given  by  Eqs.  (16)  and  (17)  of  Ref.  18. 

The  beam  current  density  i3  calculated  in  the  full  SARLAC  code  by 
accumulating  contributions  from  the  simulation  particles  onto  a  polar  mesh. 
Since  v2  ■  c,  all  simulation  particles  remain  at  constant  c  and  move  in  the 
transverse  plane  in  response  to  the  fields  calculated  from  Eqs.  (69)— (73). 
Because  Eqs.  ( 69 ) — ( 73 )  bo  not  involve  derivatives  in  z  and  because  field 


information  can  propagate  only  forward  in  5,  SARLAC  is  structured  differently 
from  conventional  simulations.  SARLAC  treats  one  c-slice  at  a  time,  thus 


reducing  the  number  of  particles  in  the  simulation  at  any  given  time  to 
-  10  .  Sach  beam  slice  is  propagated  forward  in  z  until  a  specified  maximum 
propagation  distance  z^^  is  reached;  at  this  point,  particles  are  loaded 
into  the  next  slice  at  2  -  0,  and  the  process  is  repeated. 

Nonlinear  particle  simulations  similar  to  SARLAC  have  been  developed  by 
Freeman^  and  Godfrey.^1*  The  major  novel  feature  of  the  SARLAC  model  is  its 
field  solver  which  first  does  a  predictor  step  using  the  azimuthally  averaged 
conductivity,  followed  by  two  corrector  steps.  This  process  is  fast  and  thus 
allows  a  large  number  of  Fourier  modes  to  be  retained. 

4.1  SARLAC  Particle  Simulations 

The  SARLAC  particle  code  was  designed  primarily  to  treat  the  nonlinear 
evolution  of  the  resistive  hose  instability.  Thus,  the  beam  displacements 
X(;,z)  and  Y(c,z)  in  the  full  simulation  are  calculated  self-consi3tently  from 
the  particle  dynamics.  The  results  below  are  for  a  7  kA,  5  MeV  beam  with  a 
nominal  radius  of  a0  •  1  .5  cm.  The  beam  is  injected  into  a  conducting  drift 
tube  of  radiU3  b  -  13.5  cm  filled  with  air  at  230  Torr  (p  -  0.3  atm).  The 
oeam  profile  was  tapered  at  injection  and  had  a  rise  length  of  5  -  210  cm 

(7  ns).  These  parameters  closely  resemble  the  ZTA  experimental  parameters 
reported  by  Chambers,  et  al. 

Figure  5  is  a  plot  of  the  beam  and  net  currents,  1^(0  and  I^(;,z),  at 
various  propagation  distances  z.  Initially,  In  is  everywhere  smaller  than 
Ib.  A3  the  beam  propagates  to  large  z,  the  value  of  In  rises  in  the  beam  tail 
and  soon  exceeds  .  At  z  -  270  cm  from  injection,  the  peak  net  current  is 
11,5  kA ,  as  compared  to  1^  -  6.3  kA  at  that  point  in  the  beam.  The  peak  net 
current  is  thus  70  percent  larger  than  the  beam  current. 

The  association  of  current  enhancement  with  large-amplitude  hose  motion 
is  shown  in  Figs.  7  and  3.  These  figures  plot  the  x  and  y  positions  of  a 
subset  of  simulation  particles  versus  5  at  constant  z.  As  is  weli-Khown, 1 4 


the  hoee  amplitude  increases  with  5  and  (for  a  while)  with  2.  At  z  *  1 30  cm 
(Fig.  7),  the  beam  displacements  are  of  the  order  of  the  beam  radius,  and 
current  enhancement  is  not  observed.  3y  2  -  270  cm  (Fig.  8),  the  hose  motion 
of  the  beam  is  highly  nonlinear,  and  the  beam  is  close  to  striking  the  wall. 
For  c  <  500  cm,  the  beam  radius  decreases  with  5  as  Ip  and  the  resulting  pinch 
force  build  up.  However,  in  the  tail  of  the  beam,  the  beam  expands  even 
though  Ir  exceeds  1^.  This  is  a  nonlinear  effect  which  arises  when  the  beam 
moves  out  into  regions  far  from  the  pipe  axi3  where  the  fields  generated  by 
the  front  of  the  beam  have  fallen  off. 

The  spatial  distribution  of  fields  and  currents  can  be  seen  in  contour 
plots  at  a  fixed  point  in  5  and  2.  Figure  9  is  a  contour  plot  of  the  beam 
current  density  Jb  at  5  «  630  cm  and  2  *  270  cm,  showing  that  the  beam  has 
moved  close  to  the  drift  tube  wall.  The  corresponding  plot  of  conductivity  0 
in  Fig.  10  shows  that  a  is  3pread  over  a  much  larger  area  yet  still  peaks  near 
the  pipe  axis.  Contour  plots  for  the  axial  electric  field  Ez  are  shown  in 
Fig.  11  and  display  a  characteristic  two-lobed  structure.  As  the  beam  spirals 
counter-clockwise  and  enters  virgin  air,  it  induces  a  field  which  drives 
plasma  current  opposite  to  the  beam  current.  As  the  beam  leaves  an  area,  it 
causes  the  field  to  change  sign  and  thus  drives  plasma  current  parallel  to  the 
beam  current.  The  plasma  conductivity  is  higher  in  the  latter  region. 

Current  multiplication  results  when  the  integrated  current  in  the  forward 
current  lobe  exceeds  that  in  the  return  current  lobe.  The  net  current  density 
contours,  J„  -  Jh  *  cE  ,  are  shown  in  Fig.  12.  Jn  is  spread  over  a  larger 
area  than  J>j  and  its  centroid  generally  lies  closer  to  the  pipe  axis. 

The  phenomenology  3eer.  in  this  simulation  most  closely  resembles  the 
model  in  Sec.  2.2  for  narrow  conductivity  channels.  Although  the  conductivity 
is  spread  over  a  large  area,  it  peaks  near  the  pipe  axis  and  is  low  near  the 
pipe  wall.  As  the  beam  moves  closer  to  the  wall,  the  boundary-layer  effect 
discussed  in  Sec.  3  should  emerge.  Current  gains  much  above  a  factor  of  two 
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are  unlikely,  however,  because  the  beam  is  broad  and  the  conductivity  is  low 
in  the  outer  regions.  SARLAC  runs  must  typically  be  terminated  before  the 
beam  strikes  the  wall  because  the  accuracy  of  the  field  solver  becomes  suspect 
in  this  limit. 

The  simulation  results  from  SARLAC  are  similar  to  the  ETA  experimental 
results.  In  particular,  appreciable  current  enhancement  is  generally  observed 
only  after  the  beam  has  propagated  sufficiently  far  that  large-amplitude  hose 
motion  has  occurred.  Current  enhancements  as  high  as  90  percent  (F  -  1.9) 
were  observed  on  ETA.  ^  Recent  experiments  by  Raleigh  and  Greig®  have 
confirmed  the  correlation  between  hose  displacement  and  current  enhancement 
through  use  of  a  stabilizing  soler.oidal  field.  A  solenoidal  field  of  several 
kilogauss  greatly  reduced  beam  displacement  and  eliminated  the  observed 
current  enhancement.  A  similar  effect  was  reported  by  Ipatov,  et  al.7 

4.2  Simplified  SARLAC  Simulations 

In  this  subsection  we  present  results  from  a  simplified  version  of  the 
SARLAC  code.  As  in  the  analytic  work  described  in  Secs.  2.2  and  3,  we  do  not 
use  self-consistent  beam  dynamics  but  instead  represent  the  beam  as  a  rigid 
structure  of  constant  radius  whose  vector  displacement  r&  is  specified. 

Unlike  the  analytic  work,  however,  the  simulations  do  incorporate  a  self- 
consistent  treatment  for  the  plasma  conductivity.  The  removal  of  beam 
dynamics  from  SARLAC  enables  us  to  study  the  dependence  of  current  enhancement 
on  individual  beam  and  cavity  parameters  as  distinct  from  the  dependence  on 
hose  growth  in  total. 

The  nominal  parameters  for  this  study  are:  beam  current  1^  ■  7  :<A,  beam 
duration  cmax  -  103  cm,  beam  current  rise  zr  -  210  cm,  beam  radius  a  -  1  cm, 
pipe  radius  b  -  7  cm,  and  gas  density  p  -  1  atm.  The  beam  centroid  r&  -  (X,Y) 
i3  taken  to  be  a  regular  spiral  growing  linearly  in  ;,  as  discussed  in  Ref.  14 
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for  hose  instability  in  the  beam  body:  X(c)  -  6  *  10  ’c  sin  [7.3  In  (c/3)3 

and  Y(?)  ■  6  *  10  ^5  cos  [7.3  In  (5/3)].  so  that  the  bean  displacement 

distance  r  (5)  -  6  *  10  The  current  enhancement  factor  F(;)  -  I  (;)/L  (c) 

0  n  b 

is  plotted  as  curve  (a)  in  Fig.  13.  Note  that  current  enhancement  does  not 
become  appreciable  until  approaches  half  the  wall  radius  (<;  >  500  cm). 

For  a  giver,  displacement,  the  enhancement  is  smaller  than  in  the  previous 

simulation,  due  primarily  to  a  smaller  displacement  velocity  . 

The  theoretical  basis  developed  in  Secs.  2  and  3  predicts  that  current 

enhancement  is  a  boundary  effect  arising  from  magnetic  flu/  compression 

against  either  a  metallic  pipe  or  the  "vacuum  conductivity"  radius,  whichever 

is  smaller.  Increasing  the  boundary  radius  for  a  given  beam  perturbation 

should  therefore  reduce  the  enhancement.  This  prediction  was  confirmed  by 

enlarging  the  pipe  radius  from  7  cm  to  49  cm  while  keeping  the  remaining 

simulation  parameters  unchanged.  As  seen  in  curves  (a)  and  (b)  of  Fig.  13,  F 

does  decrease  but  not  dramatically  because  the  monopole  fields  now  terminate 

at  the  vacuum  boundary,  b  *  10  cm  <  b  ,  -  49  cm. 

vac  pipe 

The  effect  of  decreasing  c.ne  beam  radius  a  from  1  cm  to  0.75  cm  i3  shown 

in  curve  (c).  The  effect  is  small,  as  expected.  The  maximum  current  gain  in 
this  case  was  F^max^  *  1.8.  The  fluctuations  in  F  seen  at  large  displacements 

in  Fig.  13  are  indications  of  field-solver  difficulty. 

An  important  effect  not  yet  discussed  is  the  polarisation  state  of  the 

hose  instability.  The  circularly  polarized  hose  mode  assumed  here  causes  a 
monotonic  increase  in  the  beam  displacement  distance  rQ(<;).  As  a  result,  the 
circuit  Inductance  L  decreases  steadily  while  the  net  current  I„  smoothly 
rises.  If  hose  growth  were  plane  polarized  (perhaps  due  to  an  asymmetric 
initiation),  the  displacement  rb,  inductance  L  and  net  current  I„  would 
oscillate,  as  seen  in  Fig.  14  where  a  beam  centroid  motion  of  X  -  0  and 
Y  ■  3  *  10  ^5  cos  [7.3  In  (?/3)]  was  imposed.  A  typical  experiment  produces 
random  hose  polarizations,  and  hence  In  should  increase  more  or  less  steadily 
but  with  some  irregular  fluctuations. 
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5. 


Conclusion 


We  have  shown  that  current  enhancement  is  a  natural  consequence  of  the 
resistive  hose  instability  for  electron  beams  propagating  in  dense  gases.  The 
degree  of  enhancement  depends  on  the  displacement  amplitude  of  the  beam 
relative  to  the  boundary  radius  where  the  monopole  electric  field  E2Q  falls  to 
zero.  The  boundary  radius  is  usually  determined  by  the  ground  returns  (e.g., 
metallic  drift  tube)  present  in  most  experiments.  If  the  return-current 
structures  are  far  from  the  beam,  current  enhancement  still  occurs  but  the 
boundary  radius  is  determined  by  the  location  where  beam-generated 
conductivity  becomes  vanishingly  small. 

The  physical  basis  for  current  enhancement  is  conservation  of  magnetic 
flux.  As  a  beam  displaces  off  axi3,  it  lowers  the  effective  inductance.  The 
circuit  responds  by  raising  the  net  current  so  as  to  conserve  magnetic  flux 
(or  energy)  within  the  conducting  plasma.  Enhancement  takes  place  provided 
the  beam  displaces  within,  a  monopole  time  -  L/R  where  L  is  the  circuit 
inductance  and  R  is  the  plasma  resistance. 

The  plasma  eddy  currents  which  produce  current  enhancement  also  drive  the 
hose  instability.  The  hose  instability  develops  on  a  plasma  dipole  decay  time 
t^.  Large  hose  growth  is  therefore  possible  in  a  monopole  time  tq  >>  td. 

Beams  undergoing  large-amplitude  hose  motion  are  thus  apt  to  produce 
substantial  current  enhancement.  Note  that  the  hose  instability,  in  contrast 
to  the  two-stream  instability,1  produces  current  enhancement  regardless  of  the 
sign  of  the  beam  charge.  A  positive-ion  beam  undergoing  hose  motion  would 
thus  induce  electric  fields  that  again  raise  the  net  current. 

Net  currents  as  large  as  three  times  the  beam  current  or  more  are 
possible  if  the  plasma  conductivity  extends  out  to  the  boundary  radius.  In 
practice,  the  current  gain  rarely  exceeds  a  factor  of  two,  due  both  to  the 


restricted  distribution  of  plasma  conductivity  and  to  the  finite  size  of  the 
beam.  Our  particle  simulation  results  showed  current  enhancements  of  -  70 
percent,  which  is  in  good  agreement  with  the  ETA  experimental  results. 
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Appendix  A:  Inductance  Formulae 


We  present  without  derivation  some  inductance  formulae  for  several 
different  conf igurations.  These  formulae  when  inserted  into  circuit  equation 
(14)  enable  one  to  compute  the  net  current  for  a  prescribed  beam  perturbation. 
The  maximum  current  enhancement  resulting  from  sudden  displacement  of  a  beam 
outside  an  on-axis  plasma  channel  is  given  as  an  example. 

The  inductance  for  an  on-axis  beam  and  plasma  of  characteristic  radius  a 
is  given  from  Eq.  (15)  by 

L  [ln( b/a)  *  a]  (Al ) 

c 

where  b  is  the  pipe  radius  and  a  is  a  parameter  that  depends  on  the  detailed 
3patial  distribution  of  the  plasma  conductivity  and  net  current.  Typically, 

a2  <  (A2 ) 


If  the  beam  and  plasma  channel  expand  self -similarly ,  the  parameter  L  in 
circuit  equation  (14)  should  be  set  to 


1  3L 

2  at’ 


(A3) 


For  a  beam  and  conductivity  channel  of  radius  a  which  are  displaced  off- 


axis  to  a  location  rw. 


W.  M .  U.  k<.L| W  -Ht  1 
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Although  the  conductivity  in  this  case  is  not  azimuthaliy  symmetric  about  the 
pipe  center  (unless  r  -  0),  Eq.  (14)  nonetheless  remains  a  reasonable 
approximation.  If  the  current  carriers  in  the  conductivity  channel  physically 
displace  with  the  beam,  one  should  set 


Otherwise,  Eq.  (A3)  should  be  used.  An  example  of  the  latter  application  is 
an  ionized  channel  consisting  of  plasma  electrons  and  ions  that  do  not  move 
laterally  but  are  created  by  the  beam  and  then  destroyed  by  recombination  as 
the  beam  moves  outward. 

If  only  the  beam  is  displaced  while  the  conductivity  channel  remains  on 
axis,  we  can  use  Eq.  (10)  to  compute  a  channel  self-inductance  given  by 


L  -  [in(b/a)  ♦  a], 
p  c2 


and  a  beam-channel  mutual  inductance  given  by 


LH  -  %  Cln(b/r. )  ♦  a'] 

b  2  b 

c 


where  in  general  a'  *  a  and  where 


r.  *  Maximum  { a ; r .  f . 

D  0 


Equation  (10)  can  then  be  rewritten  as 
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For  sudden  displacement  of  a  beam  outside  an  on-axis  plasma  channel  of 
constant  radius  a,  circuit  equation  (A9)  reduces  to  conservation  of  magnetic 
flux: 


5[L  (I  -  I.)  ♦  L.L]  -  0.  (A10 ) 

p  n  b  b  b 

If  we  ignore  the  usually  modest  changes  in  a  and  a',  relationship  (A10)  when 
coupled  to  inductance  formulae  (A6)  and  (A7)  yields 


[lr.(b/a)  +  a]  (I  -  I  )  -  tn(rK/ a)I.  . 

n  no  b  b 


(An  ) 


The  maximum  possible  current  gain  is  therefore  given  by 


F 


V1 


b 


Zr.(r.  /a) 

.  e  *  - - - 

o  Zr.(b/a)  *  a 


(A12) 


where  cQ  -  1^/1^  is  the  current  gain  prior  to  beam  displacement.  If  the 

displacement  occurs  rapidly  compared  with  a  plasma  dipole  decay  time  xd,  the 

driven  plasma  current  resides  on  the  plasma  surface  and  the  parameter  a  in  Sq. 

(A12)  should  be  set  to  zero.  For  slower  displacements,  a  is  unequal  to 

zero.  Observe  that  approximation  (A10),  and  hence  result  (A12),  is  valid  only 

if  the  displacement  occurs  rapidly  compared  with  a  monopole  decay  time, 

t  -  L  /R.  See  also  the  time-dependent  solution  given  in  Ref.  10. 

op 

Inductance  formulae  can  be  formally  given  for  broad  conductivity  channels 


but  are  generally  of  little  practical  value.  For  example,  for  uniform  a  the 
distribution  of  plasma  conductance  is  giver,  by 


Inductance  definition  (15)  then  reduces  to 


C  O  b 


~2  <’  -  ■>*> 
c 


(AH) 


where  the  mean-square  centroid  of  the  net  current  is  given  by 


,2  li 

3r- 


(A15) 


Inductance  formulae  for  the  planar  geometries  considered  in  Sec.  3  are 
similarly  straightforward.  Unfortunately,  these  expressions  depend  strongly 
on  the  spatial  distribution  of  the  electric  field  which  can  be  determined  only 
after  a  full  field  analysis.  A  circuit  analysis  is  then  unnecessary  and 


superfluous 


Appendix  S:  Solution  to  Equation  (46) 


In  this  appendix,  we  derive  the  solution  to  Eq.  (46)  via  Fourier 

analysis.  Since  the  current  enhancement  factor  F  is  independent  of  I&,  we  set 

411,/cir  -  1  and  introduce  the  dimensionless  variables  5  -  vx/Z  and 
0 

t  -  1Tt/2Tdif  to  r*ortnalize  EQa*  (46)-(48)  as 


32A. 


3S 


3A 

-  «(g  -  h(x) ) ;  0<5<ir,  0<x 


<h 


(SI  ) 


A  (0 , t)  -  A  (it, t)  -  0 
z  z 


(B2) 


A  (5,0)  -  (ir/4 )  [1  -  |  25  “  »|/ir]. 
2 


(B3) 


In  Eq.  (Bi ) , 

h(T)  -  j  *  vt  (0J*) 

where  v  i3  the  dimensionless  velocity  defined  in  Eq.  (51  )  of  the  main  text  and 
r...  is  the  characteristic  diffusion  time. 

Qlx 

Representing  the  solution  as  a  Fourier  series 

m 

A  (5,0  -  l  a  (x)  si r.  r.5,  (35) 

z  a-t 

we  note  that  the  boundary  condition  (32)  is  automatically  satisfied  and  that 


the  initial  condition  (B3)  yields 


Upon  substituting  (35)  into  (31),  we  obtain  an  equation  for  a  (t): 


2  2  rnir 

*  n  a  (t)  -  |  sin  [^L  «.  r.vx] 
ot  r.  ir  c 


(B7) 


which  may  be  easily  integrated  to  yield 


a_(t)  -  |  -r— y - =“  '  -2-J 

n  ir  _2 ,  2  ^  2» 

n  (n  *  v  ) 


{  n  sin  nh  -  r.v  cos  nh  +  e 


-n2x 


(v'Sin  nh  +  nv  cos  nh  )} 


(38) 


where  h  is  defined  in  (B-4 )  and  h  -  v/2.  The  initial  condition  (36)  has  been 

o 

used  to  obtain  (B8). 

The  current  enhancement  factor  F  i3  given  by  [cf.  (49b)]: 


F(t) 


2  l  r.a  (t) 
R’odd 


(B9) 


where  we  have  used  (B5)  to  obtain  the  last  expression.  Some  care  is  to  be 
exercised  in  evaluating  the  infinite  sum  since  the  first  term  in  the  curly 
bracket  of  (38)  contributes  to  a  nonuniformly  convergent  series  in  (39).  To 
remove  this  nonuniform  convergence,  let  us  write 


2  l  {- 

L  lir 


2  n  sin  nh 


■  4  I  (i  - 


g-Jdd  "  rX~  ‘  '  R-oio  "  n("2  *  v2) 


)  sir.  nh 


1  - 


iz 
*  l 


R’idd 


r.(n2  ♦  v2) 


sin  nh 


(BIO) 


where  we  have  used  the  identity 


4  r_,_  .  sin  3h  .  sin  5h 
_  ^la  h  *  -3  -  — 


■  j  -  1  ;  0  <  h  <  ir. 


(311  ) 


Mote  that  the  infinite  series  in  (311)  is  nonuniformly  convergent,  but 
converges  in  a  point-wise  manner  for  all  h  in  the  open  interval  (0,ir).  We  now 
substitute  (38)  into  (39),  and  use  (310)  to  write 


f(t)  .1 


-  y  ,2  2, 

it  L  ,  n(n  +  v  ) 
n”odd 


2 

{  -  v  sin  nh  -  nv  cos  nh 


-n  t  2  ,  r  nir  \ 
+  e  v  sin(^-jf 


(312) 


which  may  readily  be  compared  with  Eq.  (52)  in  the  main  text.  The  infinite 
series  in  (312)  is  absolutely  and  uniformly  convergent.  As  the  current  sheet 
is  about  to  touch  the  conducting  plate,  t  ♦  ir/2v,  h  ♦  n,  and  Eq.  (B12)  yields 


7  l  .  „2  *  v2  *  l 


(-1)  v‘e 


n-1 
n  odd 


n-1  .  , 
n  odd 


2  2 

n(n^  ♦  v  ) 


(313) 


which  is  reproduced  as  Eq.  (50)  in  the  main  text  upon  using  the  identity 


.  (  ttv  ^  ^  r  v 

tanh  I5-J  -  -  I  -5 - p  ‘ 

n-1  ..  n  ♦  v* 
n  odd 


(B14) 


This  identity  may  be  established  by  a  standard  technique  of  complex  analysis. 

Clearly,  the  method  detailed  in  this  appendix  may  readily  be  used  to 
construct  the  homogeneous  solution  Ah  and  the  driven  solution  Ad,  as  suggested 
in  Sec.  3.  For  example,  the  homogeneous  solution  A^  is  obtained  by  setting 
the  right-hand  members  of  Eqs.  (31)  and  (37)  equal  to  zero.  Then,  (35)  -  (37) 


(315) 


which  is  rewritten  as  Sq.  (53)  of  the  main  text  upon  using  (B6).  The  driven 
solution  may  be  obtained  similarly. 
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Appendix  C:  Extension  to  Two  Dimensions 


,t; 

ft 


In  this  appendix  we  demonstrate  the  insensitivity  of  the  current 

enhancement  mechanism  to  geometry  by  extending  the  analysis  of  Sec.  3  to  two 

dimensions.  We  also  demonstrate  that  our  results  are  insensitive  to  the 

initial  location  of  the  electron  beam  for  v  >>  1 . 

Consider  an  infinitesimally  thin  current-carrying  filament  ("wire") 

situated  in  a  square  channel  of  side  length  l,  inside  of  which  the  medium  i3 

characterized  by  a  constant,  uniform  electrical  oonducti vity .  We  again  assume 

that  the  wire  sits  inside  the  channel  for  a  sufficiently  long  time  prior  to  t 

-  0  that  a  steady  magnetic  field  is  already  established  at  t  »  0.  For  t  >  0, 

the  wire  is  given  a  uniform  x-ward  velocity  u. 

In  terms  of  the  normalized  variables  outlined  in  Appendix  9,  the  above 

3y3tem  is  described  by  the  vector  potential  A  (S.n.t)  whose  evolution  is 

z 

governed  by 

32A„  32A,  3A 

'  zf  '  tt  *  ir  ‘ 4(5 ' 5ln  *  V- 

95  9n 

0<5<ir,0<n<ir,T>0  (Cl) 

where  5,n  are  the  normalized  spatial  (x,y)  variables  and  t  *  1Tt/2Tdif.  as 

before.  The  initial  coordinates  of  the  wire  are  (5,n)  -  (h  ,  o  ). 

o  o 

For  r  >0,  we  assume 


h(t)  «  h  +  vt 
o 


(C2) 


and  the  instar.ter.ous  coordinates  of  the  wire  are  (£,n)  *  (h,n  ).  To 
demonstrate  the  insensitivity  to  the  geometry,  we  let  h  ,  n  be  arbitrary  in 


tills  appendix.  [Of  course,  hQ,  nQ»  h,  5,  n  are  all  restricted  to  lie  between 
0  and  it].  Equation  (Cl)  is  subject  to  the  boundary  conditions 


A  (o.n.t)  -  A  (v.n.x)  *  A U.O.x) 

Z  Z  2 


A2(5,it,T)  =*  0 


(C3) 


and  to  the  initial  condition 


AzU,n,0) 


l  l 

n*  1  aa—  1 


I_  r _ 

2  *■  2 


♦  a 


— )  sin  nho  sin 


mn  sin  n£  sin  mn 
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which  reflects  that  a  steady  magnetic  field  is  already  established  at  x  ■  0. 
The  solution  to  Eq.  (Cl )  may  be  represented  as  a  Fourier  series 


A  U.n.x)  -  l  l  a  (x)  sir.  sin  mn. 
n«i  b-1 


(C5) 


In  so  doing,  the  boundary  conditions  given  in  (C3)  are  automatically 
satisfied.  Upon  substituting  (C5)  into  (Cl),  we  obtain  the  first-order 
ordinary  differential  equation  for  a  (x): 


3amn(T)  ,2  2S  ,  ,  4 

— -  (a  +  n)a  (x)  -  -r  3in  nh(x)  sin  mn„ 

ox  an  2  0 

Tt 


ccs) 


which  may  easily  be  solved  upon  substituting  (C2)  into  (C6)  and  using  the 
initial  condition 


s  ( 0 ) 
an 


v2(n2  * 


a2) 


sin  nh  sin  an  . 
0  o 


The  current  enhancement  factor  is 


(C7) 


3  r1'  r1 

F(x)  -  1  -  yj:  J  J  d£dnAzU,n,x) 


CCS) 


which  may  also  be  expressed  as 


The  last  expression  is  analogous  to  Eq.  (39).  The  non-unifora  convergence  in 

the  infinite  series  may  be  handled  in  a  similar  manner  as  in  Eq.  (BIO). 

Equation  (C9)  may  be  used  to  calculate  the  current  enhancement  factor  as 

the  wire  is  about  to  reach  the  side  E  «  v,  [i.e.,  h  -*•  *].  This  limiting  value 

of  F  is  shown  in  Fig.  15  as  a  function  of  the  normalized  velocity  v  for 

various  combinations  of  the  initial  coordinates  (h  ,ri  ).  Mote  the 

o  o 

insensitivity  to  (h0.nQ)  and  the  similarity  between  the  curves  in  Fig.  15  and 
that  on  Fig.  1.  Thus,  the  limiting  value  of  F  *  3  is  established.  We  have 
also  shown  that  the  evolution  of  F(t),  a3  governed  by  (C9),  i3  similar  to  Fig. 
2.  Coupled  with  the  interpretation  given  in  Sec.  3.2,  we  may  conclude  that 
the  much  simpler  one-dimensional  model  provides  a  sound  physical  description 
for  current  enhancement. 


Current  gain  F  Just  before  the  current  source  touches  the 
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Figure  3: 


(b) 


Evolution  of  the  homogeneous  current  density  Jph,  in  units  of 
2 

Cirlp/Ui  ,  for  v  -  0.5  and  v  -  4.  The  arrows  at  the  bottom  mark 


Figure  4;  Evolution,  of  the  driven  current  density  Jpd,  in  units  of 

? 

cirl j^/42.  ,  for  v  -  0.5  and  v  -  4.  The  arrows  at  the  bottom  mark 
the  instantaneous  location  of  the  current  source  for  each  curve 


r.  y-.-  *ru  .npxrjVTif.  VUW 


Figure  8:  Positions  x(?)  and  y(?)  after  propagating  z  -  270  cm 


Figure  11:  Contour  plot  of  the  axial  electric  field  Ez(r,0)f  in  units  of 

mc2/eao,  at  e  =  630  cm  and  z  =*  270  cm.  The  red  and  yellow 

contours  define  regions  where  Ez  drives  return  plasma  current 

(oE  /J  <  0)  while  the  green,  blue,  and  violet  contours  are  where 
z  b 

Ez  drives  forward  plasma  current  (oEz/Jb  >  0). 
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Figure  12:  Contour  plot  of  the  r.et  current  density  J  (r,0),  in  units  of 
3  2 

me  /linage,  at  j  *  630  cm  and  z  =  270  cm.  In  almost  all  regions, 

J  /J.  >  0. 
n  b  — 


£)(// 
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and  b  <  -  »I9  cm.  For  curve  (c),  a  -  0.75  cm  and  b  .  -  7  cm. 


I 


Figure  15:  Current  gain  F  as  a  function  of  v  for  various  initial  locations 


(hQ,n0)  of  the  current  source. 
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